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SECTION  I 


INTRODUCTION 


The  release  of  certain  chemical  species  into  the  upper  atmosphere 
results  in  luminous  clouds  that  display  the  resonance  electronic-vibrational- 
rotation  spectrum  of  the  released  species.  Such  spectra  are  seen  in  rocket 
releases  of  chemicals  for  upper  atmospheric  studies  and  upon  reentry  into 
the  atmosphere  of  artifical  satellites  and  missiles.  Of  particular  interest 
in  this  connection  is  the  observed  spectra  of  certain  metallic  oxides  and 
air  diatomic  species.  From  band  intensity  distribution  of  the  spectra,  and 
knowledge  of  the  f-values  for  electronic  and  vibrational  transitions,  the 
local  conditions  of  the  atmosphere  can  be  determined  (Ref.  l). 

Present  theoretical  efforts  which  are  directed  toward  a more  complete 
and  realistic  analysis  of  the  transport  equations  governing  atmospheric 
relaxation  and  the  propagation  of  artificial  disturbances  require  detailed 
information  of  thermal  opacities  and  LWIR  absorption  in  regions  of  tempera- 
ture and  pressure  where  molecular  effects  are  important  (Refs.  2 and  3). 
Although  various  experimental  techniques  have  been  employed  for  both  atomic 
and  molecular  systems,  theoretical  studies  have  been  largely  confined  to  an 
analysis  of  the  properties  (bound-bound,  bound-free  and  free-free)  of  atomic 
systems  (Refs.  4 and  5)«  This  has  been  due  in  large  part  to  the  unavail- 
ability of  reliable  wavef unctions  for  diatomic  molecular  systems,  and  particu- 
larly for  excited  states  or  states  of  open-shell  structures.  More  recently, 
(Refs.  6-9)  reliable  theoretical  procedures  have  been  prescribed  for  such 
systems  which  have  resulted  in  the  development  of  practical  computational 
programs . 

The  theoretical  analysis  of  atmospheric  reactions  requires  the  knowl- 
edge of  the  electronic  structure  of  atoms,  ions  and  small  molecular  clusters 
of  nitrogen  and  oxygen  and  the  interaction  of  water  or  other  small  molecules 
with  these  clusters.  Knowledge  of  the  chemistry  of  metal  oxide  species, 
which  might  be  present  in  a contaminated  atmosphere,  is  also  desired. 

Because  of  the  computational  complexity  for  systems  with  large  numbers  of 
electrons,  traditional  ab  initio  theoretical  methods  are  difficult  and 
expensive  to  apply.  However,  ab  initio  computational  programs  based  on 
Slater-type  orbitals  and  incorporating  both  SCF  and  configuration  interac- 
tion analysis  have  been  developed  and  are  available  for  studies  to  chemical 
accuracy  for  systems  such  as  O2,  N2>  NO  and  their  corresponding  molecular 
ions . 
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Based  on  experience  gathered  through  detailed  calculations  for  N2j  Og, 

NO  and  many  metal  oxide  systems,  a research  program  was  undertaken  to 
theoretically  analyze  the  low-lying  excited  molecular  states  of  the  N0+  ion 
and  to  calculate  transition  probabilities  for  the  strongest  band  systems  and 
for  the  vibrational-rotational  system  of  the  ground  state  of  this  ion.  The 
calculated  data  are  presented  herein  in  both  transition  operator  and  f -number 
form  to  facilitate  their  use  in  spectroscopic  comparisons  and  as  a basis  for 
collisional  studies  of  the  N+  + 0 system. 

The  general  composition  of  this  report  is  as  follows.  In  Section  II, 
we  present  a review  of  the  current  status  of  quantum  mechanical  calculations 
for  molecular  systems.  This  is  followed  by  Section  III  which  deals  with  a 
description  of  the  mathematical  methods  which  were  employed  in  this  research. 
Included  in  Section  III  are  subsections  which  deal  with  the  construction  of 
electronic  wavefunctions , the  calculations  of  expectation  properties,  and  the 
evaluation  of  molecular  transition  probabilities.  The  calculated  results 
and  pertinent  discussions  are  presented  in  Section  IV. 
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SECTION  II 


CURRENT  STATUS  OF  QUANTUM  MECHANICAL 
METHODS  FOR  DIATOMIC  SYSTEMS 


The  application  of  quantum  mechanical  methods  to  the  prediction  of 
electronic  structure  has  yielded  much  detailed  information  about  atomic  and 
molecular  properties  (Ref.  7).  Particularly  in  the  past  few  years,  the 
availability  of  high-speed  computers  with  large  storage  capacities  has  made 
it  possible  to  examine  both  atomic  and  molecular  systems  using  an  ab  initio 
approach,  wherein  no  empirical  parameters  are  employed  (Ref.  10).  Ab  initio 
calculations  for  diatomic  molecules  employ  a Hamiltonian  based  on  the  non- 
relativistic  electrostatic  interaction  of  the  nuclei  and  electrons,  and  a 
wavefunction  formed  by  antisymmetrizing  a suitable  many-electron  function  of 
spatial  and  spin  coordinates.  For  most  applications  it  is  also  necessary 
that  the  wavefunction  represent  a particular  spin  eigenstate  and  that  it 
have  appropriate  geometrical  symmetry.  Nearly  all  the  calculations 
performed  to  date  are  based  on  the  use  of  one-electron  orbitals  and  are  of 
two  types:  Hartree-Fock  or  configuration  interaction  (Ref.  8). 

Hartree-Fock  calculations  are  based  on  a single  assignment  of  electrons 
to  spatial  orbitals,  following  which  the  spatial  orbitals  are  optimized, 
usually  subject  to  certain  restrictions.  Almost  all  Hartree-Fock  calcula- 
tions have  been  subject  to  the  assumption  that  the  diatomic  spatial 
orbitals  are  all  doubly  occupied,  as  nearly  as  possible,  and  are  all  of 
definite  geometrical  symmetry.  These  restrictions  define  the  conventional, 
or  restricted,  Hartree-Fock  (RHF)  method  (Refs.  11  and  12).  RHF  calcula- 
tions can  be  made  with  relatively  large  Slater-type  orbital  (STO)  basis  sets 
for  diatomic  molecules  with  first  or  second-row  atoms,  and  the  results  are 
convergent  in  the  sense  that  they  are  insensitive  to  basis  enlargement. 

The  RHF  model  is  adequate  to  give  a qualitatively  correct  description  of 
the  electron  interaction  in  many  systems,  and  in  favorable  cases  can  yield 
equilibrium  interatomic  separations  and  force  constants.  However,  the 
double -occupancy  restriction  makes  the  RHF  method  inappropriate  in  a number 
of  circumstances  of  practical  interest.  In  particular,  it  cannot  pro1,  ‘de 
potential  curves  for  molecules  dissociating  into  odd-electron  atoms  (e.g., 

NO  at  large  internuclear  separation),  or  into  atoms  having  less  electron 
pairing  than  the  original  molecule  [e.g.,  O2  0(^P)];  it  cannot 

handle  excited  states  having  unpaired  electrons  (e.g.,  the  state  of  0g 
responsible  for  the  Schumann-Runge  bands);  and,  in  general,  it  gives 
misleading  results  for  molecules  in  which  the  extent  of  electron  correlation 
changes  with  internuclear  separation. 

Configuration-interaction  (Cl)  methods  have  the  capability  of  avoiding 
the  limitations  of  the  RHF  calculations.  If  configurations  not  restricted 
to  doubly-occupied  orbitals  are  included,  a Cl  can,  in  principle,  converge 
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to  an  exact  wave function  for  the  customary  Hamiltonian.  However,  many  Cl 
calculations  have  in  fact  been  based  on  a restriction  to  doubly-occupied 
orbitals  and  therefore  retain  many  of  the  disadvantages  of  the  RHF  method 
(Ref.  8).  The  use  of  general  Cl  formulations  involves  three  considerations, 
all  of  which  have  been  satisfactorily  investigated:  the  choice  of  basis 

orbitals,  the  choice  of  configurations  (sets  of  orbital  assignments),  and 
the  specific  calculations  needed  to  make  wavefunctions  describing  pure  spin 
states  (Ref.  6).  The  first  consideration  is  the  art  associated  with  quantum 
mechanical  electronic  structure  calcinations.  Many  methods  (iterative  NSO, 
perturbation  selection,  first  order  Cl,  etc.)  have  been  advocated  for  the 
optimum  choice  of  configurations . There  are  no  firm  rules  at  present  and 
the  optimum  choice  is  a strong  function  of  the  insight  of  the  particular 
research  investigator.  The  last  consideration,  proper  spin  and  symmetry 
projection,  has  proved  difficult  to  implement,  but  computer  programs  have 
( been  developed  for  linear  projection  algebra  at  this  Center,  and  the  Cl 
method  has  been  found  of  demonstrable  value  in  handling  excited  states  and 
dissociation  processes  which  cannot  be  treated  with  RHF  techniques. 

Either  of  the  above  described  methods  for  ab  initio  calculations 
reduces  in  practice  to  a series  of  steps,  the  most  important  of  which  are 
the  evaluation  of  molecular  integrals,  the  construction  of  matrix  elements 
of  the  Hamiltonian,  and  the  optimization  of  molecular  orbitals  (RHF)  or 
configuration  coefficients  (Cl).  For  diatomic  molecules,  these  steps  are  all 
comparable  in  their  computing  time,  so  that  a point  has  been  reached  where 
there  is  no  longer  any  one  bottleneck  determining  computation  speed.  In 
short,  the  integral  evaluation  involves  the  use  of  ellipsoidal  coordinates 
and  the  introduction  of  the  Neumann  expansion  for  the  interelectronic 
repulsion  potential  (Ref. 13);  the  matrix  element  construction  depends  upon 
an  analysis  of  the  algebra  of  spin  eigenfunctions  (Ref.  1^);  and  the  orbital 
or  configuration  optimization  can  be  carried  out  by  eigenvalue  techniques 
(Refs.  15,  16).  All  the  steps  have  by  now  become  relatively  standard,  and 
can  be  performed  efficiently  on  a computer  having  32,000  to  65,000  words 
of  core  storage,  a cycle  time  in  the  microsecond  range,  and  several  hundred 
thousand  words  of  peripheral  storage. 

Both  the  RHF  and  Cl  methods  yield  electronic  wavefunctions  and  energies 
as  a function  of  the  internuclear  separation,  the  RHF  method  for  one  state, 
and  the  Cl  method  for  all  states  considered.  The  electronic  energies  can 
be  regarded  as  potential  curves,  from  which  may  be  deduced  equilibrium 
internuclear  separations,  dissociation  energies,  and  constants  describing 
vibrational  and  rotational  motion  (including  anharmonic  and  rotation- 
vibration  effects).  It  is  also  possible  to  solve  the  Schrodinger  equation 
for  the  motion  of  the  nuclei  subject  to  the  potential  curves,  to  obtain 
vibrational  wavefunctions  for  use  in  transition  probability  calculations . 

The  electronic  wavefunctions  themselves  can  be  used  to  estimate  dipole 
moments  of  individual  electronic  states,  transition  moments  between 
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different  electronic  states,  and  other  properties.  While  all  of  the 
calculations  described  in  this  paragraph  have  been  carried  out  on  some 
systems,  the  unavailability  of  good  electronic  wavefunctions  and  potential 
curves  has  limited  actual  studies  of  most  of  these  properties  to  a very 
small  number  of  molecules. 

A few  studies  illustrating  the  scope  of  the  current  work  in  this  field 
are  cited.  Exhaustive  RHF  calculations  have  been  reported  for  first  and 
second-row  hydrides  (Refs.  IT  and  l8)  for  most  first-row  diatomic  molecules 
(Ref.  19),  and  for  many  other  molecules  containing  second-row  atoms  (Ref. 
20).  Configuration-interaction  calculations  restricted  to  double  occupancy 
are  illustrated  by  the  work  of  Das  and  Wahl  on  Lig  and  F2  (Refs.  8 and  21). 
They  determined  an  optimum  choice  of  orbitals  for  a small  number  of  config- 
urations designed  to  permit  proper  description  of  the  dissociation  products, 
and  obtained  highly  satisfactory  potential  curves.  Davidson  (Refs.  22  and 
23)  has  carried  out  doubly  occupied  Cl  studies  with  very  large  numbers  of 
configurations  to  gain  more  insight  into  the  description  of  correlation 
energy.  The  largest  of  the  Cl  studies  not  restricted  to  doubly  occupied 
orbitals  have  been  carried  out  at  UTRC.  The  62  states  of  O2  dissociating 
into  low-lying  atomic  oxygen  states  have  been  described  in  a qualitatively 
consistent  manner  (Ref.  2U).  Many  of  these  states  involve  several  impaired 
electrons,  and  the  success  of  the  treatment  depended  critically  upon  the 
inclusion  of  all  types  of  pertinent  configurations  and  upon  proper  handling 
of  the  spin.  Similar  work  on  the  102  low-lying  states  has  now  been 
completed  (Ref.  25).  In  striking  contrast  to  the  recent  progress  in 
obtaining  electronic  energies  and  wavefunctions,  rather  few  calculations 
of  electronic  transition  moments  have  been  attempted.  Among  the  few  studies 
in  this  area  is  the  work  of  Michels  on  He  H+  (Ref.  26)  and  of  Henneker  and 
Popkie  (Ref.  27)  on  diatomic  hydrides  using  Hartree-Fock  wavefunctions. 

More  recently  a theoretical  program  for  calculating  band  strengths  for  the 
systemsNg  (1  PS),  0g  (SR),  and  NO  (p)  has  been  carried  out.  The  results 
of  this  program  indicate  that  reliable  band  strengths  (10  to  25  percent) 
can  be  calculated,  provided  a Cl  approach  is  employed  (Ref.  28).  Similar 
work  has  been  performed  for  band-to-bai  d transition  probabilities  for  the 
systems  LiO,  A10,  FeO,  UO,  U0+  and  TiO  (Ref . 29  ). 
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SECTION  III 


METHOD  OF  APPROACH 


A.  Electronic  Structure 

A spin-free  nonrelativistic  Hamiltonian  is  employed  in  the  Bom-Oppenheimer 
approximation.  In  systems  containing  atoms  as  heavy  as  N or  0,  this  approxi- 
mation is  quite  good  for  low-lying  molecular  states.  For  a diatomic  molecule 
containing  n electrons,  the  approximation  leads  to  an  electronic  Hamiltonian 
depending  parametrically  on  the  internuclear  separation  R: 


Za  ZB 


I 


(1) 


1=  I 


l=l 


'A  in  iB 


t >)>  I 


where  Z.  and  Z are  the  charges  of  nuclei  A and  B,  and  ri.  is  the  separation 
a ^ a A . 

of  electron  i,  and  nucleus  A.  The  Hamiltonian  H is  in  atomic  units  (energy  in 

Hartrees,  length  in  Bohrs). 

Electronic  wavefunctions , i|/(r)  , are  made  to  he  optimum  approximations 
to  solutions  of  the  Schrodinger  equation 


^(R)t(R)=  E (R)*(R) 

by  invoking  the  variational  principle 

/V(RU(R)V(R)tfr 


(2) 


(3) 


8w(R)  = S 


jV(R)'i'(R)a'T 


The  integrations  in  Eq.  (3)  are  over  all  e1ectronic  coordinates  and  the 
stationary  values  of  W(R)  are  approximations  to  the  energies  of  states  des- 
cribed by  the  corresponding  i|/(r).  States  of  a particular  symmetry  are  studied 
by  restricting  the  electronic  wavefunction  to  be  a projection  of  the  appro- 
priate ancmlar  momentum  and  spin  operators.  Excited  electronic  states  corres- 
ponding to  a particular  symmetry  are  handled  by  construction  of  configuration 
interaction  wavefunctions  of  appropriate  size  and  form. 
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The  specific  form  for  if/( R)  may  be  written 


*(R)  = Z c„*u.  (R) 

H- 

where  each  i//  (R)  is  referred  to  as  a configuration,  and  has  the  general 

r 

structure 


(M 


V<R>  = ^S  n ( r.i,R )eH 


i^i 


(5) 


where  each  if/  . is  a spatial  orbital,  /7 is  the  antisymmetrizing  operator, 

Q g is  the  spin-projection  operator  for  spin  quantum  number  S,  and  0^  is 
a product  of  a and  /3  one-electron  spin  functions  of  magnetic  quantum  number 
M.  No  requirement  is  imposed  as  to  the  double  occupancy  of  the  spatial  orbi- 
tal, so  Eqs.  (M  and  (5)  can  describe  a completely  general  wavefunction. 


In  Hartree-Fock  calculations  if/  is  restricted  to  a single  if/^  which  is 
assumed  to  consist  as  nearly  as  possible  of  doubly-occupied  orbitals.  The 
orbitals  if/^  are  then  selected  to  be  the  linear  combinations  of  basis  orbi- 
tals best  satisfying  Eq.  (3).  Writing 


= 1 

1/ 


Xv 


(6) 


The  a^.  are  determined  by  solving  the  matrix  Hartree-Fock  equations 


\ Fx,a,,  = €i?  sx,a„  <each  M 


(7) 


where  e . is  the  orbital  energy  of  if/  ■ 

3.  " r*  1 

The  Fock  operator  F has  been  throughly  discussed  in  the  literature 
(Ref.  31)  and  depends  upon  one-  and  two-electrcn  molecular  integrals  and  upon 
the  avi.  This  makes  Eq.  (7)  nonlinear  and  it  is  therefore  solved  iteratively. 
UTRC  has  developed  programs  for  solving  Eq.  (7 ) for  both  closed  and  open- 
shell  systems,  using  basis  sets  consisting  of  Slater-type  atomic  orbitals. 
Examples  of  their  use  are  in  the  literature  (Ref.  6). 

In  configuration  interaction  calculations,  the  summation  in  Eq.  (U  ) 
has  more  than  one  term,  and  the  c^  are  determined  by  imposing  Eq.  (3),  to 
obtain  the  secular  equation 


X ( HfjLv 


ws 


/xv 


) Cv  = 0 


(eoch  /x ) 


(8) 


7 


where 


=/'f'/!( R)«^  (R)  ^1/  (R)  dr 
S/ii/  (R)'I'i/  (R)  dr 


Equation  (7)  is  solved  by  matrix  diagonalization  using  either  a modified 
Givens  method  (Ref.  15)  or  a method  due  to  Shavitt  (Ref.  16). 

The  matrix  elements  H|jV  and  Sqv  may  be  reduced  by  appropriate  operator 
algebra  to  the  forms 

€P<8m  #sp|0^>^n  ^i(xi-R)  J (r)p  n "*'^i (i j, r))>  do) 

V=£  €P  <0M|0S  p|e^><(^n  (Xj.R)  P n (£1  ,r^>  (n) 

where  P is  a permutation  and  its  parity.  The  sum  is  over  all  permutations. 

^ | @ 5 P | 9m>  is  a "Sanibel  coefficient"  and  the  remaining  factors  are 

spatial  integrals  which  can  be  factored  into  one-and  two-electron  integrals. 

If  the  are  orthonormal,  Eqs.  (io)  and  (ll)  become  more  tractable  and  the 
H^v  and  S|jVmay  be  evaluated  by  explicit  methods  given  in  the  literature 
(Ref.  14).  Computer  programs  have  been  developed  for  carrying  out  this  pro- 
cedure, and  they  have  been  used  for  problems  containing  up  to  Uo  total  elec- 
trons, 10  unpaired  electrons,  and  several  thousand  configurations. 

The  Cl  studies  described  above  can  be  carried  out  for  any  orthonormal 
set  of  ’Vi  for  which  the  molecular  integrals  can  be  calculated.  Programs 
developed  by  UTRC  made  specific  provision  for  the  choice  of  the  </^  as 
Slater-type  atomic  orbitals,  as  symmetry  molecular  orbitals,  as  Hartree- 
Fock  orbitals,  or  as  more  arbitrary  combinations  of  atomic  orbitals. 

In  summary,  computer  programs  have  been  developed  which  are  capable  of 
carrying  out  all  the  steps  needed  to  make  diatomic  Hartree-Fock  or  Cl  studies 
for  closed-  or  open-shell  systems,  including  excited  states,  based  on  Slater- 
type  orbitals.  These  studies  lead  to  electronic  energies  and  wavefunctions 
as  a function  of  the  intemuclear  separation. 
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B. 


Vibrational  and  Rotational  Properties 


For  an  electronic  state  described  by  if)  ( R)  and  W(r),  the  relative  motion 
of  the  nuclei  is,  in  the  Bom-Oppenheimer  approximation,  subject  to  a poten- 
tial having  at  internuclear  separation  R the  value  W(R).  By  considering  the 
quantum  mechanics  of  the  nuclear  motion,  it  is  possible  from  W(R)  to  calculate 
vibrational  and  rotational  energy  levels.  It  is  convenient  to  report  the 
vibration-rotation  structure  in  terms  of  parameters  r , Dg,  &>eXg»  Be>  ae>  etc., 
which  are  also  available  by  standard  reduction  of  experimental  data.  Such  an 
analysis  can  be  routinely  carried  out  for  bound  electronic  states,  using  a 
Dunham  analysis  computer  program  which  has  been  deposited  in  the  Quantum 
Chemistry  Program  Exchange  (QCPE  113)* 

From  W(R)  it  is  also  possible  to  obtain  vibrational  wavefunctions  by 
numerical  integration  of  the  radial  Schrodinger  equation  for  the  nuclear  motion. 
A computer  program  which  carried  out  this  calculation  by  the  Numerov  proce- 
dure has  been  developed.  Any  W(R)  can  be  handled  since  the  program  fits  it  by 
a spline  technique.  This  program  has  been  used  for  several  years  on  a variety 
of  problems;  a typical  application  has  been  to  excitation  transfer  in  colli- 
sions of  normal  and  metastable  He  atoms  (Ref.  30  ).  The  input  W(R)  can  be 
the  calculated  potential  corresponding  to  the  electronic  wavefunctions  (R)  or 
it  can  be  derived  using  an  RKR  procedure  (Refs.  31  and  32).  The  input  can  be 
either  the  basic  B(v)  and  G(v)  data  or  the  derived  experimental  spectroscopic 
constants  based  on  this  data.  Programs  for  implementing  the  RKR  procedure 
have  been  described  in  the  literature  (Ref. 33). 
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C.  Transition  Probabilities 


The  electronic  and  vibrational -rotational  wavefunctions  of  a pair  of 
states  can  be  used  to  calculate  transition  probabilities.  If  two  molecular 
states  are  separated  in  energy  by  an  amount  AE^  = hcv  (h  = Planck's  constant, 
c = velocity  of  light,  v = frequency  in  wave  numbers),  the  semiclassical  theory 
of  radiation  (Refs.  34  and35)yields  for  the  probability  of  a spontaneous 
transition  from  an  upper  state  n to  a lower  state  m 


A nm  = x 


AE  nm 

h4c3 


Jnm 

Qn 


(12) 


Here  A is  the  Einstein  coefficient  for  spontaneous  transition  from  level 
nm 

n-*m,  g is  the  total  degeneracy  factor  for  the  upper  state 
n 


9n  = (2-  8cr,A,)(2S,+  l)(2J'  + l) 


(13) 


and  S is  the  total  strength  of  a component  line  in  a specific  state  of 
polarization  and  propagated  in  a fixed  direction.  A related  quantity  is 
the  mean  radiative  lifetime  of  state  n defined  by 

= Z A nm  (l4  ) 

n men 

the  summation  being  over  all  lower  levels  which  offer  allowed  connections. 

The  intensity  of  the  emitted  radiation  is 


Inm  = AEnmNnAnm 


(15) 


where  N is  the  number  density  in  the  upper  state  n.  This  analysis  assumes 
that  all  degenerate  states  at  the  same  level  n are  equally  populated,  which 
will  be  true  for  isotropic  excitation.  The  total  line  strength  can  be 
written  as  the  square  of  the  transition  moment  summed  over  all  degenerate 
components  of  the  molecular  states  n and  m: 


S 


nm 


Mji 

• .)  1 


(16) 


10 


where  j and  i refer  to  all  quantum  numbers  associated  collectively  with  upper 
and  lower  electronic  states,  respectively. 

In  the  Born-Oppenheimer  approximation,  assuming  the  separability  of 
electronic  and  nuclear  motion,  the  wavefunction  for  a diatomic  molecule  can  be 
written  as 

*VJMA=  ^el  (£•  R)  '/'v  (R)  V'jmA  X-  4> ) (17) 

where  ijj1  (r,  R)  is  an  electronic  wavefunction  for  state  i at  fixed  inter- 
nuclear6separation  R,  tfj  (R)  is  a vibrational  wavefunction  for  level  v and 
if,  (0,  X ,<£)  refers  to  the  rotational  state  specified  by  electronic  angular 
momentum  A,  total  angular  momentum  J and  magnetic  quantum  number  M.  The  rep- 
resentation is  in  a coordinate  system  related  to  a space-fixed  system  by  the 
Eulerian  angles  ( 6 , X,<f>  ).  The  transition  moment  M _ can  be  written,  using 
the  wavefunction  given  by  Eq.  ( 17) , as  ^ 


M. 


= Mj.A'«'{Se  + Kn}+ivA-«"  dT'dl«dr- 


(18) 


The  subscripts  e,  v and  r refer  to  the  electronic,  vibrational  and  rotational 
wavefunctions  and  and  JMp  are  the  electronic  and  nuclear  electric  dipole 
moments,  respectively.  Integration  over  the  electronic  wavefunction,  in  the 
Born-Oppenheimer  approximation,  causes  the  contribution  of  the  nuclear  moment 
Mn  to  vanish  for  : = .1 . The  electronic  dipole  moment  can  be  written  (Refs. 

35  and  36) in  the  form 

Me  = -Z  { y (®,x  ,+)  (19) 


where  the  primed  coordinates  refer  to  the  space  fixed  system,  the  coordinates 
r refer  to  a molecule-fixed  system  and.2>(0, X,*/*)  is  o.  group  rotation  ten- 
sor whose  elements  are  the  direction  cosines  related  to  the  Eulerian  rotation 
angles  ($,X,<f>).  Using  bracket  notation,  Eqs.  (l8)  and  (19)  can  be  combined 
to  yield  for  the  transition  moment 


jv'j'A'*' 

M j j = M iv"j"A,'m" 


'£e~ 


k 


<J'A'M'  5 (0,x,*)|  J*A"M y 


(20) 
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The  matrix  elementsCJ'  A'  M'  |2>(0,  X ,<f>)  | J"  A " M">determine  the  group  selec- 
tion rules  for  an  allowed  transition  and  have  been  evaluated  for  many  types  of 
transitions  (Refs.  37,  36  and  39)  Summing  Eq.  (20)  over  the  degenerate  mag- 
netic quantum  numbers  M'  and  M"  we  have  from  Eq.  (l6) 


_ _nv'j'A'  _ ^J'A'  nv' 
>>nm-  imv"j"A"  ~ ^J"A"  Pmv" 


j f A * 

where  £ is  the  Honl -London  factor  (Refs.  4o  and  1*1 ) and 

J"  A " 


<v'Ne4v> 


is  the  band  strength  for  the  transition.  Combining  Eqs.  (l3),  (15)  and  (2l), 
we  have  for  the  intensity  of  a single  emitting  line  from  upper  level  6: 


f ac’wV  ]4cnvVA' 
nv'j'  _ A ^mv*j"A" 


Inm=  1 mv«j»  = 3 Nj/ 


Ti4C3wn(2j'+|) 


where  N " xs  the  number  density  in  the  upper  rotational  state  J'  and  to  = 

(2-  (2S'  + l)  is  the  electronic  degeneracy.  Taking  an  average  vaSue  of 

En  v,',  J,',  for  the  whole  hand,  Eq.  (23)  can  be  summed  to  yield  the  total  inten- 

IT  V J 

sity  in  the  (v',v")  band: 


,nv'  nv'j'  . 4 , [aeRJ  Pmy» 

Imv  3 Nv  r.4c3u*n 


where  N ' = ji  Nji  is  the  total  number  density  in  the  upper  vibrational  level 
v*  and  #here  we  make  use  of  the  group  summation  property 


z iw'-Hi 


1 


Comparing  Eq.  (15)  and  (24),  we  have  for  the  Einstein  spontaneous  transition 
coefficient  of  the  hand  (v',v") 


Tac  nv<  l3r>nV,„ 

_ _4_  [AE  mV»J  PmVw 
3 k W.i. 


nv1 

*mv 


(26) 


tl*C  3W| 


Similarly,  the  lifetime  of  an  upper  vibrational  level  v'  of  state  n can  be 
written 


. nv 


* -L  V™" 


(27  ) 


where  the  summation  runs  over  all  v"  for  each  lower  state  m.  Equation  (26  ) 
can  be  cast  in  the  computational  form 

»SC'.(s.c-).  U'4'ff.»JSl>  [ae^.,  (ou)]  Vi"  too) 


(28  ) 


where  £E^  y"  and  pn  are  in  atomic  units.  It  is  also  often  convenient  to 
relate  the  transition  probability  to  the  number  of  dispersion  electrons  needed 
to  explain  the  emission  strength  classically.  This  number,  the  f-number  or 
oscillator  strength  for  emission,  is  given  by 


fnm,W  = 


me3  h2 


2e 


:Ue"Vv'"] 


A nv 
*mv 


(29) 


The  inverse  process  of  absorption  is  related  to  the  above  development 
through  the  Einstein  B coefficient.  Corresponding  to  Eq.  (15),  we  have  for  a 
single  line  in  absorption 


'mn  r 
? Ax  J 


v line  (vVjV) 


K (>/)  d 1/  - hj/mnNmBmn 


(30) 


where  K(v)  is  the  absorption  coefficient  of  a beam  of  photons  of  frequency 
and 

nv'j'A' 

nv'J'A'  2 n 

Bmn=  Bmv"j"A"  = 3f)*c  wm(2J"+l) 


(31) 
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is  the  Einstein  absorption  coefficient  for  a single  line.  Summing  over  all 
lines  in  the  band  (v",  v’),  assuming  an  average  band  frequency,  we  obtain 


= NV«  PmV'1 

V 3Tl2Ccum  ^mV 


where  N " = N „ is  the  total  number  density  in  the  lower  vibrational  state 
v".  Corresponding  to  Eqs.  (28)  and  (29)  we  can  define  an  f-number  or  oscil- 
lator strength  for  absorption  as 


f mn.v'V  = 


— nv' 

2m  AEmV»  nv« 

3T.V«.m  PmV' 


In  computational  form,  Eq.  (33)  becomes 

f f_  2 Ae|),Vv«(q.u.)  nv ' i3k 

' mn,v  v * 3 PmV"  u) 

where  AEn  v,’,  and  pn  V„  are  in  atomic  units.  Combining  Eqs.  (26)  and  (29)  and 
comparing  with  Eq.10^),  we  see  that  the  absorption  and  emission  f-numbers  are 
related  by 


fmn,v"v'  = (S^)  fnm.v'  v" 


Some  caution  must  be  observed  in  the  use  of  f-numbers  given  either  by  Eq.  (29) 
or  ( 33)  since  both  band  f-numbers  and  system  f-numbers  are  defined  in  the  lit- 
erature. The  -onfusion  arises  from  the  several  possible  band  averaging  schemes 
that  can  be  identified. 

An  integrated  absorption  coefficient  (density  corrected)  can  be  defined 
from  Eq.  (32)  as 


sv'v'  = ^T  Cv,,=  Nv"Bv':v'  ('-exP 


■hc^y*y'  \ hl/yW  yf 


(36) 


where  the  exponential  factor  corrects  for  stimulated  emission.  Equation  (36) 
can  be  written  in  terms  of  the  absorption  f-number  as 


e ne 

Svys 


exp 


- hci'y* 
kT 


j fmn,v" 


v' 


(37) 


Using  hvc/k  = 1.1+3880  cm-K°,  we  obtain  a computational  formula  for  the 
integrated  absorption  coefficient  as  S , (cm-^  . atm  = 


v V’ 

Sy'v'  (cm-2  atm-l)  = 

„ /ov*  Ift\  / Nu*\  / 143880  1/V"v' (cm-1)  \ z 

*■ 3795  * 107  TO)  (U7  ) ('-exp T I' 

The  total  integrated  absolution  is  found  from 


(38) 


v' 


ST0TAL=£  £ SV*,V' 

V"  v' 


(39  ) 


where,  under  normal  temperature  conditions,  only  the  first  few  fundamentals 
and  overtones  contribute  to  the  summations. 


The  developments  given  above  are  rigorous  for  band  systems  where  an 
average  band  frequency  can  be  meaningfully  defined.  Further  approximations, 
however,  are  often  made.  For  example,  the  electronic  component  of  the  dipole 
transition  moment  can  be  defined  as 


fcO  ) 


This  quantity  is  often  a slowly  varying  function  of  R and  an  average  value  can 
sometimes  be  chosen.  Equation  (22  ) can  then  be  written  approximately  in  fac- 
tored form  as 

Pmv"  ~ Qv'v*  £ I ^ ji  (R)|  O+l) 

where  q , „,  the  square  of  the  vibrational  overlap  integral,  is  called  the 

Fran  ck-£onc[on  factor.  Q is  evaluated  at  some  mean  value  of  the  inter- 
nuclear  separation  R.  ^addition , it  is  sometimes  possible  to  account  for  a 
weak  R-dependence  in  Mf  by  a Taylor  series  expansion  of  this  quantity  about 
some  reference  value , usually  referred  to  the  (0,0)  band.  We  have 
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1 


- fiyP [ l + a(R-Ra/3)  + b(R-Ra/9)2+  ] 


0*2  ) 


Substituting  into  Eq.  (22)  and  integrating  yields 


Pmv*  ~ P v'v"  X | ^ji  [ l + a (Rv'v"  ~ Ra/3)  + b(Rv'v"  - Ra/3)2+  ] | 0*3) 


where 


<V'|(R- RQi3)|  V'^> 


0*1* ) 


is  the  R-centroid  for  the  transition  and 


, _ ,2  <y'|  (R-Ra>3)2lv,5> 

(Rv'v"-Ra/?)  = v>> 


(1*5) 


is  the  R -centroid.  Note  that  this  last  \ erm  differs  (to  second  order)  from 
the  square  of  the  R-centroid.  An  alternate  procedure  can  be  developed  by 
evaluating  Eq.  (1*0 ) at  each  R-centroid,  Rv,  v«.  Then 


Pmv"  85  Qv'v"£  |/?ji(Rv'v")| 
1 » J 


(1*6) 


Equation  (1*6)  assumes  that  the  vibrational  wavefunction  product  ib  \l/  behaves 
like  a delta  function  upon  integration. 


V v 


Vv  Vv,= S(R~  Rv'v") 


(1*7 ) 
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The  range  of  validity  of  Eq.  ( l^g)  is  therefore  questionable,  particularly  for 
band  systems  with  bad  overlap  conditions  such  as  oxygen  Schumann-Runge . The 
range  of  validity  of  the  R-centroid  approximation  has  been  examined  by  Frazer 
(Ref.  1+2). 


The  final  step  in  calculating  transition  probabilities  is  the  determination 
of  /?  (R),  the  electronic  dipole  transition  moment,  for  the  entire  range  of 

internuclear  separations,  R,  reached  in  the  vibrational  levels  to  be  consid- 
ered. This  can  be  expressed  in  terms  of  the  expansion  of  Eq.  ( 1+ ) as 


^ji(R)=Z  cu#  cj,  <>yi(R)|  Me|  ^„(R)>  (l+8) 

/i.1/  ^ 

where  cj^  and  c*  are  coefficients  for  if/ ^ and  , respectively. 

An  analysis  similar  to  that  yielding  Eqs.  (lO)  and  (ll)  gives 

I «p  <^ms|  0sP|^<^n  ^k(ik.R)|Mep|  n^K^K.R)) 

(1+9) 

The  spatial  integral  in  Eq.  (1+9)  reduces  to  one-electron  integrals  equivalent 
to  overlap  integrals,  and  the  evaluation  of  Eq.  (1+9)  can  be  carried  out  by  the 
same  computer  programs  used  for  Eq.  (ll).  Programs  for  evaluating  /?  (R)  in 

Eq.  (1+8)  have  been  developed  at  UTRC  and  examples  of  their  application 
have  appeared  in  the  literature  (Ref.  26). 

For  perturbed  electronic  systems,  the  transition  dipole  moment  will  have 
a strong  R-dependence  and  R-centroid  or  other  approximations  will  be  invalid. 

A direct  evaluation  of  Eq.  (22)  would  therefore  be  required  using  the  fully- 
coupled  system  of  electronic  and  vibrational  wavefunctions  to  properly  account 
for  the  source  of  the  band  perturbations. 
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SECTION  IV 


DISCUSSION  OF  RESULTS 


The  absorption  coefficient  for  the  N0+  fundamental  vibration-rotation 
band  has  recently  been  measured  by  Bien,  et  al.  (Ref.  43).  The  reported 
integrated  value  is  l40  cm-2  atm--*-  in  conflict  with  previous  experimental 
estimates  of  500  cm-2  atm--*-  (Ref.  1+1)  and  a recent  theoretical  value  of  89 
cm-2  atm--*-  (Ref.  H5).  Since  the  dominant  atmospheric  molecular  ion  above 
100  km  is  NO  , its  radiation  characteristics  and  interactions  with  electrons 
and  other  atomic  and  molecular  species  are  of  primary  importance  in  under- 
scanding  the  detailed  chemistry  of  the  atmosphere.  Relatively  little  is 
known  about  the  excited  electronic  states  of  the  N0+  ion.  Accurate  potential 
curves  are  available  for  the  ground  X"*£+  state  and  for  theA-4l  and  a 3j’+ 
excited  states  of  N0+  (Ref.  16).  Gilmore  also  gives  some  experimental  esti- 
mates of  other  excited  states  of  N0+.  More  recently  Thulstrup,  et  al. 

(Refs.  17  and  18)  have  assigned  additional  excited  states  of  N0+  based  on 
photoelectron  spectrum.  A summary  of  these  data  is  shown  in  Fig.  1.  A 
correlation  diagram  for  some  of  the  low-lying  valence  excited  states  of 
N0+  is  shown  in  Table  I.  In  addition  to  these  valence  excited  states 
there  exists  a large  manifold  of  Rydberg  states  to  the  N+(3p)  + 0+(2+S) 
limit.  The  energies  of  the  dissociation  limits  of  the  low-lying  elec- 
tronic states  of  N0+  are  shown  in  Table  2.  It  is  obvious  that  only  a 
small  fraction  of  the  low-lying  excited  states  for  this  system  have 
been  experimentally  investigated. 


Few  theoretical  studies  have  been  carried  out  for  the  excited  states 
of  this  system.  Lefebvre-Brion  and  Moser  (Refs.  19  and  50 ) have  studied 
the  lowest-lying  % and  3n  states  in  an  SCF  framework.  More  recently, 
Thulstrup  and  Ohrn  (Refs  47  and  48)  have  examined  low-lying  singlet  and  trip- 
let states  of  E+,  £-,  II  ana  a symmetry.  Their  studies  excluded  higher  spin 
and  angular  momentum  states  which  may  be  important  in  perturbation  analysis. 
Also,  their  calculations  for  ^n,  -*^  and  \ symmetries  were  truncated  in 
expansion  length  owing  to  their  computer  program  limitations.  Their  pre- 
dicted locations  for  these  states  are  not  in  good  agreement  with  the  best 
experimental  estimates. 


A complete  and  quantitatively  uniform  ab  initio  Cl  study  of  all  the  low- 
lying  valence  states  of  N0+  was  carried  out  as  part  of  the  research  conducted 
under  this  program.  This  study  was  similar  to  that  carried  out  for  the  102 
low-lying  valence  states  of  the  nitrogen  molecule  (Ref.  25).  Many  of  the 
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states  are  similar  since  N0+  is  isoelectronic  with  Ng-  Empirically  adjusted 
potential  energy  curves,  based  on  both  the  ab  initio  calculations  and  the 
available  experimental  spectroscopic  data,  have  been  derived  and  studies  of 
the  radiation  characteristics  of  the  strongest  electronic  transitions  have 
been  carried  out. 

Figures  2-17  illustrates  the  low-lying  valence  state  potential  energy 
curves  for  the  N0+  system.  These  curves  have  been  calculated  using  a full 
Cl  within  a minimum  STO  basis.  The  configuration  sizes  and  number  of  states 
calculated  are  shown  in  Table  3 . 

These  calculations  have  the  property  of  correct  dissociation  to  atomic 
limits  and,  as  previously  found  for  the  N2  system,  a gradual  and  uniformly 
increasing  correlation  error  with  decreasing  internuclear  separation.  This 
error  arises  from  inadequacy  of  the  basis  set  to  fully  describe  polarization 
effects  between  the  core  and  valence  electrorls.  This  error  is  quantitatively 
similar  for  the  excited  low-lying  states  of  a given  symmetry  and  our  calculated 
potential  curves  have  therefore  been  empirically  adjusted  to  yield  agreement 
with  RKR  curves  constructed  from  the  limited  known  specroscopic  data.  This 
empirical  adjustment  takes  the  form: 

&E(R)  - E(R)calc.  “ E^R^exp. 

We  find  that  is  a strongly  increasing  function  of  R for  short  internu- 
clear separations  but  is  nearly  independent  of  the  spin  state  and  only 
weakly  dependent  of  the  angular  momentum  of  the  excited  states.  Thus  a 
nearly  universal  correction  curve  (error  ~ .2  eV)  exists  for  most  of  the  N0+ 
system  of  valence  excited  states. 

Our  summary  of  known  experimental  spectroscopic  data  for  the  low-lying 
bound  valence  states  of  N0+  is  given  in  Table  4.  The  corresponding  potential 
curves  are  given  in  Fig.  1.  The  most  interesting  curves  from  the  standpoint 
of  the  charge  transfer  reaction  are  those  of  and  symmetry.  An  avoided 
crossing  between  ^T.+  I and  II  is  apparent  and  this  could  be  the  most 
important  charge  transfer  channel.  Another  observation  is  that  the  5p  I 
state  connecting  with  N+  + 0 crosses  the  '>E+  I connecting  to  N + 0+  at  or 
near  the  dissociation  limit  of  N+  + 0.  This  means  that  a large  low-energy 
(<  .5  eV)  cross-section  may  be  found  for  this  reaction.  Further  studies  of 
these  symmetries  are  indicated. 


k 
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We  have  carried  out  an  analysis  of  the  vibrational-rotational  radiation/ 
absorption  characteristics  for  the  ground  X z!+  state  of  N0+.  Calculated 
oscillator  strengths  for  the  vibrational- rotational  system  are  given  in 
Table  5.  The  corresponding  calculated  Einstein  coefficients  are  given  in 
Table  6.  These  are  based  on  the  calculated  dipole  moment  function.  This 
function  is  similar  to  that  given  by  Billingsley  (Ref.  1+5)  but  exhibits  a 
somewhat  steeper  slope  near  the  equilibrium  separation  of  the  ion.  The 
long  range  behavior  of  the  dipole  moment  function  for  N + 0 can  be  written 
as 


where 


therefore, 


P-(R)  - + g “ 16  Z c.m. 

R -*  » 


c.m. 


= [i  - 


mn 


Mn  + M 


] K 


0 


n(R)  - -.0311933  R 

R -•  oo 


Here  the  dipole  moment  is  defined  relative  to  a center-of-mass  coordinate 
system  and  we  note  that  jj,  (R ) is  not  gauge -invariant  for  an  ionized  system. 


This  dipole  moment  functions  yields  somewhat  higher  integrated  band 
absorption  coefficients  than  those  reported  by  Billingsley.  Our  calculated 
value  for  the  fundamental  band  is  168  cm-2  atm-1  which  can  be  compared  with 
an  experimental  lower  bound  (>  l40  cm-2  atm'1)  as  reported  by  Bien  et  al. 
(Ref.  1+3).  A complete  summary  of  the  calculated  integrated  absorption 
coefficients  is  given  in  Table  7. 

The  N0++  system  is  illustrated  in  Fig.  18  using  the  data  of  Thulstrup, 
et  al.  (Ref.  1+8).  A correlation  diagram  for  the  low-lying  valence  excited 
states  of  N0++  is  shown  in  Table  8.  The  energies  of  the  dissociation 
limits  of  the  low-lying  electronic  states  of  N0+*+are  shown  in  Table  9. 

The  configuration  sizes  of  the  symmetries  for  NO  are  presented  in  Table  10. 
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These  studies  of  the  electronic  structure  of  N0+  have  included  calcula- 
tions of  the  potential  energy  curves  to  large  internuclear  separations  where 
the  adiabatic  atomic  connections  can  be  made.  These  curves  thus  constitute 
a necessary  data  base  for  describing  the  charge  exchange  reaction: 

N+  + 0 -*  N + 0+ 

This  reaction  is  of  significant  concern  in  the  overall  DNA  chemistry  pro- 
gram of  atmospheric  modeling.  The  low  energy  (<  5 eV)  cross-section  for 
this  reaction  is  extremely  difficult  to  obtain  experimentally  and  the  pre- 
sent studies  should  establish  the  background  information  needed  for  a theo- 
retical description  of  this  reaction.  Such  a program  would  be  a natural 
follow-on  to  the  present  research  program. 
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TABLE  1 


LOW-LYING  MOLECULAR  STATES  OF  NO  AND  THEIR 
DISSOCIATION  LIMITS 


Dissociation  Limit 
(Calculated  Order) 


Molecular  States 


N(\)  + 0+(4Su) 

V(D,  V(l),  V(l), 

N+(3Pg)  + 0(3Pg) 

V(2),  V(2),  V(2), 
V(D,  V(l),  V(l), 

3n(2),  5n(2), 
^(1),  3A (1 ) , 5A(1) 

n(2du)  + o+(4su) 

V(i),  V(i),  3n(i), 

5n(i),  3a (i ) , 5a(i) 

N(4SU)  0+(2Du) 

V(1),  V(i),  3n  (i ) , 

5n(i),  3a (i ) , 5a(i) 

n(2pu)  + 0+(4su) 

V(i),  V(i),  \{i), 

5n(i) 

N+(1Dg)  + 0(3Pg) 

V(i),  V(2),  3n(3), 

3A(2),  3f(l) 

n(4su)  + o+(2pu) 

V(i),  V(d, 

3n(i),  5n(i) 

TABLE  1 (Cont'd) 


LOW-LYING  MOLECULAR  STATES  OF  NO+  AND  THEIR 
DISSOCIATION  LIMITS 


Dissociation  Limit 
(Calculated  Order) 

N+(3P  ) + 0(!d  ) 

g 8 


N+(1S  ) + 0(3P  ) 

g 8 

N(\)  + 0+(\) 


N+(1D  ) + 0(1D  ) 

g g 


n+(3p  ) + 0(1S  ) 

g g 


Molecular  States 


V(1), 

V(2), 

3n(3), 

3A(2), 

3$(1) 

V(i), 

3n(i) 

V(3), 

V(3), 

V(2) 

V(2), 

^(4), 

3n(4), 

1A(3), 

3A(3),  1 

*(2), 

3*(2),  ' 

Lr(i),  3 

r(i) 

V(3), 

V(2), 

1n(4), 

1A(3),  ' 

L$(2),  1 

r(i) 

V(l), 

3n(i) 
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TABLE  2 


ENERGIES  OF  NITROGEN  AND  OXYGEN  ATOMIC  AND  IONIC  STATES  REPRESENTING 
DISSOCIATION  LIMITS  OF  LOW-LYING  NO+  STATES 


Atomic  and 
Ionic  States 
(Experimental 
order ) 

Total  Energy  (Hartrees) 

Energy  Relative 
to  4SU  + 4SU  (eV) 

Calc. 

Exp. 

n(4su)  + o+(4su) 

-128.42277 

0.000 

0.000 

N+(3Pg)  + 0(3Pg) 

-128.31777 

2.857 

0.916 

n(2du)  + o+(4su) 

-128.30317 

3.255 

2.383 

N+(\)  + 0(3Pg) 

-128.23804 

5-027 

2.815 

N+(3Pg)  + O^Dg) 

-128.22386 

5.413 

2.883 

n(4su)  + o+(2du) 

-128.28190 

3.833 

3.325 

n(2pu)  + o+(4su) 

-128.26192 

4.377 

3.575 

N+(1Dg)  + O^g) 

-128.14412 

7.583 

4.782 

N+(1Sg)  + 0(3pg) 

-128.17437 

6.759 

4.968 

n(4su)  + o+(2?u) 

-128.22932 

5.264 

5.016 

N+(3Pg)  + OC^-Sg) 

-128.14286 

7.617 

5.105 

n(2du)  + o+(2du) 

-128.16229 

7.088 

5.708 

U7 


TABLE  3 


CONFIGURATION  SIZES  AND  NUMBER  OF  STATES  FOR 

VCI 

CALCULATIONS  OF  VARIOUS 

SYMMETRIES  FOR  NO 

Symmetry 

No.  Cfgs 

No.  of  States 

(full  Cl) 

(to  and  including 
N+3P  + 0 XS 
d i s s §c  iat  ion  ®1  im  it 

296 

9 (+),  5 (-) 

-E 

392 

10  (+),  11  (-) 

5l 

126 

5 (+),  3 (-) 

?£ 

10 

1 ( + ) 

in 

240 

10 

3n 

328 

18 

5n 

96 

7 

7n 

8 

!a 

142 

7 

3a 

170 

10 

5a 

43 

3 

^A 

1 

_* ** 

48 

4 

3$ 

56 

4 

5$ 

8 

_*** 

1r 

10 

2 

3r 

6 

1 

* First  state  dissociates  to  a high  limit  N+(5su)  + 0(3p  ) 

**  First  state  dissociates  to  a high  limit  N(^P  ) + 0+(1+P  ) 

***  First  state  dissociates  to  a high  limit  N+(3§  ) + 0 (3p  ) 

u 6 


TABLE  h 


Calculated  using  experimental  constants  of  Thulstrup,  et  al.  (Ref.  1+8). 


N0+  (X  12+  - X 1Z'+) 


v’/v" 

n 

0 

1 

2 

3 

4 

5 

u 

1 

7.072-06 

2 

6.163-08 

1.408-05 

3 

1.936-09 

1.828-07 

2.100-05 

4 

1.293-10 

6.714-09 

3.686-07 

2.775-05 

5 

1.295-12 

7.394-10 

1.455-08 

6.136-07 

3.435-05 

6 

2.333-13 

2.388-11 

1.535-09 

2.636-08 

9.151-07 

4.078-05 

7 

1.109-12 

1.410-11 

8.597-11 

2.548-09 

4.254-08 

1.272-06 

8 

7.310-13 

9.871-12 

3.187-11 

1.601-10 

3.888-09 

6.401-08 

9 

8.058-12 

2.088-12 

6.4i8-li 

4.181-11 

2.709-10 

6.195-09 

10 

5.324-12 

3.967-12 

1.900-11 

1.304-10 

1.415-11 

3.694-10 

v’/v" 

6 

7 

8 

9 

10 

0 

1 

2 

3 

It 

5 

6 

7 U. 703-05 

8 1.691-06  5.305-05 

9 8.944-08  2.176-06  5.884-05 

10  9.491-09  1.191-07  2.729-06  6.435-05 
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TABLE  6 


Calculated  Einstein  Coefficients  (kv<v«  , 
for  the  X - X System  of  N0+ 


v'/v" 

0 

1 

2 

3 

4 

5 

0 

1 

2.592+01 

2 

8.912-01 

5.019+01 

3 

6.211-02 

2.569+00 

7.273+01 

4 

7.272-03 

2.093-01 

5.034+00 

9.340+01 

5 

1.121-04 

4.039-02 

4.406-01 

8.140+00 

1.123+02 

6 

2.868-05 

2.009-03 

8.146-02 

7.754-01 

1.179+01 

1.293+02 

7 

1.829-04 

1.683-03 

7.023-03 

1.313-01 

1.214+00 

1.590+01 

8 

1.552-04 

1.581-03 

3.695-03 

1.270-02 

1.944-01 

1.773+00 

9 

2.133-03 

4.303-04 

9.975-03 

4.703-03 

2.084-02 

3.004-01 

10 

1.714-03 

1.019-03 

3.800-03 

1.966-02 

1.544-03 

2.756-02 

v'/v" 

6 

7 

8 

9 

10 

0 

1 

2 

3 

It 

5 

6 

7 1.447+02 

8 2.050+01  1.583+02 

9 2.402+00  2.558+01  1.702+02 

10  4.463-01  3.100+00  3.108+01  1.804+02 


TABLE  7 

Total  Integrated  Absorption  Coefficient 
for  N0+  (X  V) 


-2 

Absorption  Coefficient,  cm  atm 


Temperature 

First 

Second 

Third 

(°K) 

Fundamental 

Overtone 

Overtone 

Overtone 

Total 

100. 

168.278 

1.467 

0.046 

0.003 

169.794 

273.16 

168.278 

1.467 

0.046 

0.003 

169.794 

300. 

168.278 

1.467 

0.046 

0.003 

169.794 

500. 

168.276 

1.470 

0.046 

0.003 

169.795 

1000. 

168.223 

1.569 

0.050 

0.004 

169. 846 

1500. 

168.087 

1.811 

0.06l 

0.005 

169.964 

2000. 

167.886 

2.133 

0.076 

0.007 

170.102 

2500. 

167.631 

2.501 

0.095 

0.009 

170.236 

3000. 

167.328 

2.894 

0.118 

0.012 

170.352 

4000. 

166.590 

3.730 

0.172 

0.013 

170.510 

5000. 

165.682 

4.6o6 

0.236 

0.026 

170.550 

8000. 

161.637 

7.260 

0.46l 

0.057 

169.415 

10000. 

157.481 

8.776 

o.6o4 

0.077 

166.938 
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TABLE  8 


T— 


++ 

Low-Lying  Molecular  States  of  NO  and  Their 
Dissociation  Limits 


Dissociation  Limit 
(Experimental  Order) 

N+(3P  ) + o+(4s  ) 
g u 

n++(2p  ) + 0(3P  ) 
u g 

N+( “4)  ) + 0+(4S  ) 
g u 

n+(3p  ) + o+(2d  ) 
g u' 

N++(2P  ) + 0(4)  ) 

u g 

N+(1S  ) + 0+(^S  ) 
g u' 

n+(3p  ) + o+(2p  ) 

g U 

N+(1D  ) + 0+(2D  ) 
g u 

n++(2p  ) + o(1s  ) 
u g 

n+(5su)  + o+(4su) 

n+(4> ) + o+(2p  ) 

g u 

n(4s  ) + 0++(3P  ) 
u g 


Molecular  States 


Vd),  V<d,  Vd), 

2n(i),  4n(i),  6n(i) 

2z'(2),  V(2),  2E+(i),  4 Z+(  1 ) 

2n(2),  4n(2),  2a(i),  4a(i) 

V(i),  4n(i),  4a(i) 

2E‘(1),  V(l),  V(2),  V(2), 

2n(3),  4n(3),  2a(2),  4a(2), 
2*U),  hHD 

V(i),  2 E+( 2 ) , 2n(3),  2a(2), 

2*(1), 


2I"(2),  4z"(2),  V(l),  V(l), 

2n(2),  4n(2),  2a(i),  4a(i) 
V(3),  V(2),  2n(4),  2a(3), 

2H2),  2r(i) 

V(i),  2n(i) 

2I+(l),  V(i),  CV(1),  8E+(1) 

V(i),  2 z+  ( 2 ) . 2n(3),  2 a( 2 ) , 

2*(1) 

V(i),  V(i),  V(D,  2n(i), 

4n(i),  6n(i) 


TABLE  9 


Energies  of  Nitrogen  and  Oxygen  Ionic  States  Representing 
Dissociation  Limits  of  Low-Lying  NO  States 


Ionic  States 
(Experimental  Order) 

N (JP  ) + 0 ( S ) 
g u 

++.2  . ,3  . 

N ( P ) + o(JP  ) 
u g 

N+(1D  ) + 0+(4S  ) 
g u 

n+(3p  ) + o+(2d  ) 

g U 

N++(2P  ) + oC1!)  ) 
u g 

N+(1S  ) + 0+(1+S  ) 
g u 

N+(3P  ) + o+(2p  ) 
g u 

) + o+(2d  ) 

g u 

n++(2p  ) + o(1S  ) 
u g 

n(^s  ) + 0 ( S ) 
u u 

N+( “4)  ) + 0+(2P  ) 

g u 

N(  S ) + 0 ( P ) 

u g 


Total  Energy  (Hartrees) 

-127.93207 

(-127.20712)* 

-127.85233 

-127.79139 

(-127.1132)* 

-127.76201 

-127.69747 

-127.70185 

(-127.03220)* 

-127.84272 

-127.60793 

(-127.09326)* 


Energy  Relative 

to  N+(V  ) + oYs  ) 
g(eV) 


Calc. 

0.000 

0.000 

(19.727)* 

1.452 

2.170 

1.899 

3.828 

3.325 

(22.283)* 

3.M9 

4.628 

4.052 

6.384 

5.016 

6.265 

5.223 

(24.487)* 

5.641 

2.431 

5.847 

8.820 

7-015 

(22.825)* 

7.002 

* 

non-opt  imized  calculation 
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Concord  Sciences 
P.  0.  Box  119 
Concord,  MA  01742 
Att:  E.  A.  Sutton 

Kent  State  University 
Department  of  Chemistry 
Kent , OH  44242 

Attn:  W.  C.  Fernelius 

Ecometrics,  Inc. 

718  Willow  Glen  Road 
Santa  Barbara,  CA  93105 
Attn:  N.  G.  Utterback 

Epsilon  Laboratories,  Inc. 

4 Preston  Court 
Bedford,  MA  01730 
Attn:  J.  Dalchines 

C.  Accardo 
H.  Miranda 


ESL,  Inc. 

495  Java  Drive 
Sunnyvale,  CA  94086 
Attn:  J.  Marshall 

W.  R.  Bell 


General  Electric  Company 
Space  Division 
Valley  Forge  Space  Center 
Goddard  Blvd  King  of  Prussia 
P.0.  Box  8555 
Philadelphia,  PA  19101 
Attn:  R.  H.  Edsall 

P.  Zavitsanos 
J . Burns 
T . Baurer 
F.  Alyea 

M.  H.  Bortner  Space  Sci  Lab 
M.  Linevsky 


General  Electric  Company 
Temp-Center  for  Advanced  Studies 
8l6  State  Street  (P.0.  Drawer  QQ) 
Santa  Barbara,  CA  93102 
Attn:  B.  Gambill 

W.  S.  Knapp 
DASIAC  A.  Feryok 
DASIAC 
T.  Stephens 
D.  Chandler 

General  Motors  Corporation 
Delco  Electronics  Division 
Santa  Barbara  Operations 
6767  Hollister  Avenue 
Goleta,  CA  93017 
Attn:  R.  I.  Primich 

General  Research  Corporation 
P.0.  Box  3587 
Santa  Barbara,  CA  93105 
Attn:  H.  H.  Lewis 

J.  Ise  Jr. 

General  Research  Corporation 
Washington  Operations 
Westgate  Research  Park 
7655  Old  Springhouse  Road 
Mclean,  VA  22101 

Attn:  T.  M.  Zakrzewski 

Geophysical  Institute 
University  of  Alaska 
Fairbanks,  Alaska  99701 
Attn:  B.  J.  Watkins 

T.  N.  Davis 
H.  Cole 
N . Brown 

R.  Parthasarathy 
D.  J.  Henderson 

J.  S.  Wagner,  Physics  Dept 

Govt  Fiscal  Relations  & Patent  Off 
275  Admin  Building  AG-70 
University  of  Washington 
Seattle  WA  98l95> 

Attn:  R.  Geballe 

K.  C.  Clark 


63 


Department  of  Defense  Contractors 
(Cont'd) 

Harvard  University- 
Department  of  Atmospheric  Sciences 
Pierce  Hall 
Cambridge  MA  02138 
Attn:  M.  B.  McELroy 

HSS , Inc. 

2 Alfred  Circle 
Bedford,  MA  01730 
Attn:  D.  Hansen 

M.  P.  Shuler 
H.  Stewart 

Institute  for  Defense  Analyses 
400  Army-Navy  Drive 
Arlington,  VA  22202 
Attn:  H.  Wolfhard 

E.  Bauer 
J.  Bengston 
P.  A.  Selwyn 

IRT  Corporation 
P.  0.  Box  81087 
San  Diego,  CA  92111 
Attn:  J.  A.  Rutherford 

R.  H.  Neynaber 
D . A . Vroom 

Johns  Hopkins  University 
Applied  Physics  Laboratory 
Johns  Hopkins  Road 
Laurel  MD  20810 

Attn:  Document  Librarian 


Kaman  Sciences  Corporation 
P.0.  Box  7H63 

Colorado  Springs,  CO  80933 
Attn:  B.  J.  Bittner 

P.  Jesson 

F.  H.  Shelton 

IIT  Research  Institute 
10  West  35th  Street 
Chicago,  111.  60616 

Attn:  Technical  Library 


Lockheed  Missiles  and  Space  Company 
3251  Hanover  Street 
Palo  Alto,  CA  94304 

Attn:  M.  Walt  Dept  52-10 

J.  B.  Regan  d/52-12 
B.  M.  McCormac  Dept  52 -l4 
James 

L.  Kulander  Dept  52-14 
J . Kumer 

R.  D.  Sears  Dept  52-l4 

G.  Johnson  Dept  52-12 
C.  Gunton  Dept  52-l4 


T. 

J. 


R. 

R. 


Maryland,  University  of 
Administration  Building 
College  Park,  MD  21030 
Attn:  T.  J.  Rosenberg 

R.  J.  Munn  Dept  of  Chem 
J.  Vanderslice  Chem  Dept 
Dept  of  Physics  H.  Griem 

Maxwell  Laboratories,  Inc. 

9244  Balboa  Avenue 
San  Diego,  CA  92123 
Attn:  V.  Fargo 

McDonnell  Dougles  Corporation 
5301  Bolsa  Avenue 
Huntington  Beach,  CA  92647 
Attn:  J.  Moule 

W.  Olsen 
A.  D.  Goedeke 

Photometries,  Inc. 

442  Marrett  Road 
Lexington,  MA  02173 
Attn:  I.  L.  Kofsky 

Physical  Dynamics  Inc. 

P.0.  Box  1069 
Berkeley,  CA  94701 
Attn:  A.  Thompson 

J.  B.  Workman 
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Department  of  Defense  Contractors 
(Cont'd) 


Physical  Sciences,  Inc. 
30  Commerce  Way 
Woburn,  MA  01801 


Attn:  R.  L.  Taylor 

K.  Wray 


Physics  International  Company 

2700  Merced  Street 

San  Leandro,  CA  9^577 

Attn:  DOC  CON  FOR  Tech  Lib 


Pittsburgh,  University  of 
Of  the  Comwlth  Sys  of  Higher  Educ 
Cathedral  of  Learning 
Pittsburgh,  PA  15213 
Attn:  M.  A.  Biondi 

E.  Gerjudy 

F.  Kaufman 
W.  L.  Fite 


Princeton  Univ.  The  Trustees  of 
Forres tal  Campus  Library 
Box  36 

Princeton  University 
Princeton,  NJ  08540 
Attn:  A.  J.  Kelly 


R&D  Associates 
P.0.  Box  9695 

Marina  Del  Rey,  CA  90291 
Attn:  B.  Gabbard 

R.  Latter 
F.  Gilmore 

C.  H.  Humphrey 

D.  Dee 

R.  P.  Turco 
H.  A.  Dry 
R.  G.  Lindgren 
R.  E.  Lelevier 
A.  L.  Latter 


R&D  Associates 

1815  N.  Ft.  Myer  Drive 

11th  Floor 

Arlington,  VA  22209 
Attn:  H.  J.  Mitchell 

J.  M.  Rosengren 


Rand  Corporation,  The 
1700  Main  Street 
Santa  Monica,  CA  90406 
Attn:  C.  Crain 


Rice  University 
Department  of  Space  Science 
Houston,  TX  77058 
Attn:  R.  F.  Stebbings 


Midwest  Research  Institute 
425  Volker  Boulevard 
Kansas  City,  MO  64llO 

Attn:  Physics  Sec  T.  H.  Milne 
Physics  Sec  F.  T.  Greene 


Minnesota,  University  of 
2030  University  Avenue,  S.E. 
Minneapolis,  MN  554l4 
(No  Classified  to  this  Address) 
Attn:  J.  R.  Winkler 


Minnesota,  University  of 
Morris  Campus 
Morris,  MN  56267 
Attn:  M.  Hersch 


Mission  Research  Corporation 
735  State  Street 
Santa  Barbara,  CA  93101 
Attn:  R.  Hendrick 

C.  L.  Longmire 
R.  Kilb 

D.  Archer 

M.  Scheibe 
R.  N.  Wilson 
D.  Sowle 

D.  Sappenfield 
W.  F.  Crevier 

N.  J.  Carron 
R.  Bogus ch 
P.  Fischer 


Department  of  Defense  Contractors 
(Cont'd) 


Mitre  Corporation,  The 
P.0.  Box  208 
Bedford,  MA  01730 
Attn:  S.  A.  Morin  M/s 

J.  N.  Freedman 

National  Academy  of  Sciences 
Attn:  Nat’l  Materials  Advisory  Board 
2101  Constitution  Avenue 
Washington,  D.C.  204l8 
Attn:  W.  C.  Bartley 

J.  R.  Sievers 
E.  R.  Dyer  (2  copies) 

New  York  at  Buffalo,  State  Univ.  of 
Buffalo,  NY  14214 
Attn:  G.  0 Brink 

Pacific-Sierra  Research  Corp. 

1456  Cloverfield  Blvd 
Santa  Monica,  CA  90404 
Attn:  E.  C.  Field  Jr 


Riverside  Research  Institute 
80  West  End  Avenue 
New  York,  NY  10023 
Attn:  G.  Glaser 

J.  B.  Minkoff 
D.  H.  Koppel 

Science  Applications,  Inc. 

P.0.  Box  2351 
La  Jolla,  CA  92038 
Attn:  D.  A.  Hamlin 

B.  Myers 
R.  W.  Lowen 

Science  Applications,  Inc. 
Huntsville  Division 
2109  W.  Clinton  Avenue 
Suite  700 

Huntsville,  AL  35805 
Attn:  N.  R.  Byrn 

Professor  Chalmers  F.  Sechrist 
155  Electrical  Engineering  Bldg. 
University  of  Illinois 
Urbana , IL  6l801 


Pennsylvania  State  University 
Industrial  Security  Office 
Room  5 5 Old  Main  Bldg 
University  Park,  PA  16802 
Attn:  L.  Hale 

J.  S.  Nisbet 


Rice  University,  William  Marsh 
P.0.  Box  2692 
Houston,  TX  77001 
Attn:  Industrial  Sec  Super 

J.  Chamberlain 


Stanford  Research  Institute 
333  RavenswooG  Avenue 
Menlo  Park,  CA  94025 

Attn:  W.  G.  Chestnut  (2  copies) 

B.  R.  Gasten 
A.  L.  Whitson 
A.  M.  Peterson 
R.  White 
R.  L.  Leadabrand 
J.  T.  Moseley 
T.  G.  Slanger 
D.  J.  Hildenbrand 
L.  L.  Cobb 
J.  R.  Peterson 
J.  Lomax 
R.  D.  Hake  Jr. 


G.  Black 
F.  T.  Smith 
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Department  of  Defense  Contractors 
(Cont'd) 


Stanford  Research  Institute 
306  Wynn  Drive,  N.W. 

Huntsville,  AL  35805 
Attn:  M.  Morgan 

Systems,  Science  and  Software,  Inc. 
P.0.  Box  1620 
La  Jolla,  CA  92038 
Attn:  R.  C.  Sklarew 


Technology  International  Corporation 
75  Wiggins  Avenue 
Bedford,  MA  01730 
Attn:  W.  P.  Boquist 

TRW  Systems  Group 
One  Space  Park 
Redondo  Beach,  CA  90278 
Attn:  r,  k.  Plehuch 

J.  F.  Friichtenicht/Rl-1196 
Tech  Inf or  Center/S-1930 

R.  Watson  Rl/1096 

United  Technologies  Research  Center 
East  Hartford,  CT  06l08 

Attn:  B.  H.  Bullis  & H.  H.  Michels 

Utah  Stat  _ University 
Logan,  UT  84321 

Attn:  W.  M.  More  Chem  Dept 

D.  Burt 

C.  Wyatt 

D.  Baker 

K.  Baker 


Visidyne,  Inc. 

19  Third  Avenue 
North  West  Industrial  lark 
Burlington,  MA  01803 
Attn:  R.  Waltz 

J.  W.  Carpenter 
T.  C.  Degges 
0.  Manley 


Wayne  State  University 
1064  Mackenzie  Hall 
Detroit,  MI  48202 

Attn:  P.  K.  Rol  (Chem.  Eng.  & Mat.Sci.) 

R.  H.  Kummler 

Wayne  State  University 
Dept,  of  Physics 
Detroit,  MI  48202 
Attn:  W.  E.  Kauppila 

Yale  University 
New  Haven,  CT  06520 

Attn:  Engrg  Dept.  G.  J.  Schultz 

Continued  from  Department  of  the  Army 

Commander 

Harry  Diamond  Laboratories 
2800  Powder  Mill  Road 
Adel phi , MD  20783 

Attn:  M.  H.  Weiner  DRXDO-TI 

: DRXDO-NP  (2  copies) 

Commander 
Trasana 

White  Sands  Missile  Range  NM  88002 
Attn:  R.  E.  Dekinder,  Jr. 

Di rector 

U.  S.  Army  Ballistic  Res.  Labs. 

Aberdeen  Proving  Ground,  MD  21005 
Attn:  M.  D.  Kregel 

DRXBR-AM  J.  S.  Batten 
J.  C.  Mester 
Tech  Lib  E.  Baicy 
J.  A.  Vanderhoff 
DRXBR-AM  D.  Snider 
DRXBR-CA  T.  J.  O'Brien 
J.  Heimerl  DRXRD  -BSP 
I.  L.  Chidsey 
DRXBR-AM  G.  E.  Keller 


Chief  C-E  Services  Division 
TJ.  S.  Army  Communications  CMD 
Pentagon  Rm  1B269 
Washington,  D.  C.  20310 
Attn:  CC-0PS-CE 


Commander 

US  Army  Electronics  Command 
Fort  Monmouth,  N.J.  07703 

Attn:  DRSEL  -PL-ENV  H.  A.  Bomke 


Department  of  Defense  Contractors 
(Cont ' d) 

California  Institute  of  Technology- 
Jet  Propulsion  Laboratory 
U800  Oak  Park  Grove 
Pasadena,  CA  91103 
Attn:  J.  AJello 

California,  State  of 
Air  Resources  Board 
9528  Telstar  Avenue 
El  Monte,  CA  91731 
Attn:  L.  Zafonte 

Denver,  University  of 
Colorado  Seminary 
Denver  Research  Institute 
P.  0.  Box  10127 
Denver,  CO  80210 

Attn:  Sec.  Officer  for  Mr.  Van  Zyl 

Sec.  Officer  for  D.  Mur Cray 


Howard  University 
Department  of  Chemistry 
Washington,  D.  C.  20059 
Attn:  W.  Jackson 

JAYC0R 

lUOl  Camino  Del  Mar 
Del  Mar,  CA  920ll+ 

Attn:  H.  B.  Levine 

Johns  Hookins  University 
3^th  & Charles  Streets 
Baltimore,  MD  21218 
Attn:  J.  J.  Kauftnan 

Lowell,  University  of 
Center  for  Atmospheric  Research 
U50  Aiken  Street 
Lowell,  Mass.  Ol85*+ 

Attn:  G.  T.  Best 

Southern  California,  Univ.  of 
University  Park 
Los  Angeles,  CA  90007 
Attn:  S.  W.  Benson 

Stanford  Research  Institute 
l6ll  North  Rent  Street 
Arlington,  VA  22209 
Attn:  W.  W.  Beming 
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